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Abstract 

We study a solvable QCD-like toy theory, a generalization of the Nambu- 
Jona-Lasinio model, which implements chiral symmetries of light quarks 
and heavy quark symmetry. The chiral symmetric and chiral broken 
phases can be dynamically tuned. This implies a parity doubled heavy- 
light meson system, corresponding to a (0~, 1~) multiplet and a (0 + , 1 + ) 
heavy spin multiplet. Consequently the mass difference of the two mul- 
tiplets is given by a Goldberger-Treiman relation and qa is found to be 
small. The Isgur-Wise function, £(u>), the decay constant, and other 
observables are studied. 
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I. Introduction 

In recent years, QCD applied to systems containing a single very massive quark, where 
one can imagine the limit M — > oo to be a reasonable physical approximation, has 
been the subject of considerable attention [1-5]. The pseudoscalar and vector mesons 
containing one very massive and one light quark become degenerate in the M — > oo 
limit, due to a heavy quark spin symmetry again valid to 1/M. Moreover, Isgur and 
Wise [1] pointed out that transition amplitudes, such as weak decays, involving heavy 
quarks are described by a flavor independent function of the invariant difference in 
4-velocities, £(V • v), and therefore a heavy quark spin-flavor symmetry, SU(2Nf) 
exists, valid to order 1/M. Georgi has given a useful field theoretic construction of 
this limit [5] , and has studied the consequences and phenomenological applications of 
the theory, such as the computation of the QCD anomalous dimension which controls 
the perturbative evolution of £(V • v) with scale for v' • v < 1. 

For many purposes one must also implement the chiral symmetries of the light 
quarks, in addition to the heavy quark symmetry. The heavy quark (HQ) and chiral 
light quark (LQ) symmetries together control the interactions of heavy-light (HL) 
mesons with pions and TT-mesons, etc. Several authors have written down model 
independent chiral Lagrangians which involve these symmetries at the meson level 
[6-12]. A number of studies of the phenomenological applications of these chiral 
Lagrangians have been undertaken, such as the computation of the chiral log radiative 
corrections to £(V • v), [7] associated with ££7(3) x £77(3) breaking terms, the study 
of radiative and meson decays of heavy mesons [11], and chiral dynamics including 
the effects of excited heavy mesons [12]. 

The chiral Lagrangian introduced by Wise [6] represents a straightforward imple- 
mentation of the heavy quark and light flavor symmetries in the nonlinear current 
form. In this form one need only identify the linear flavor symmetries, like isospin or 
£77(3), and the chiral effective Lagrangian, to leading order in the momentum expan- 
sion, is automatically determined, up to an unspecified axial vector coupling constant 
qa- This effective Lagrangian is then manifestly invariant under the usual global flavor 
symmetries, and the full set of chiral transformations are local gauge-like transfor- 
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mations which are functionals of the pions. The underlying chiral representations of 
the heavy mesons need never be specified. Model independent approaches are clearly 
the most reliable way in which to minimally implement the physical symmetries. 

We may wish, however, to go closer to the underlying chiral dynamics than the 
model independent approaches allow. We may pose additional questions within dy- 
namical models which can reveal additional physical consequences to the real world. 
For example, is there a more primitive chiral form of the Lagrangian in which the ex- 
plict chiral representations of the heavy mesons are identified? A related question in 
the broken phase is: what is the analogue of the Goldberger-Treiman relation in the 
heavy meson system, i.e., what receives mass from the chiral condensate's mass gap? 
In the case of the nucleon-meson system we can similarly write the chiral Lagrangian 
in the nonlinear current form, never having to specify the precise chiral representa- 
tions of nucleons. However, if we ask for the linear chiral form we also know the 
answer: the left-handed (right-handed) nucleon is assigned to a (0, |), ((^,0)) repre- 
sentation under SU(2) R x SU(2) L . Most of the nucleon mass arises from the chiral 
condensate, or the VEV of £ which is (|, \). We know this because the Goldberger- 
Treiman relation yields the pion-nucleon coupling constant in terms of the nucleon 
mass Qnntt ~ ™>n/ fn, and Ga ~ 1- 

In the case of heavy mesons, however, it is clear that the meson mass arises 
primarily from the mass of the heavy constituent quark, such as the 6-quark, and the 
chiral mass gap is a perturbation. Our question then is related to the outcome of a 
gedanken experiment: what happens to the heavy-light meson spectrum if we could 
somehow restore the chiral symmetry, maintaining the other features of confining 
QCD? While the nucleon mass goes to zero in this gedanken limit, leaving degenerate 
(approximately) massless left- and right-handed states, the heavy meson masses must 
remain (approximately) unaffected. Yet, the explicit linear chiral symmetry SU (2)^ x 
577(2) # must somehow be realized in the heavy meson mass spectrum in this limit. 
This leads to the conclusion that the ground-state must become doubly degenerate 
with even and odd parity mesons B 1 and B 2 respectively, and these must form the 
SU (2) l x SU (2)r representations in the linear combinations (Bi±B 2 )/ y/2. Therefore, 
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the breaking of the chiral symmetry leads to a mass gap between these parity partners 
and associated pionic transitions between parity partners will occur (at the end of 
Appendix B we give a schematic discussion of parity doubling). 

It is difficult to imagine that a simple potential model can capture this phe- 
nomenon. The chiral symmetry limit is relativistic, and the chiral symmetry breaking 
is a dynamical rearrangement of the vacuum. Thus, the naive picture of a heavy me- 
son as a boundstate of a heavy quark and a constituent quark will miss those aspects 
of the physics which involve the necessary mixing of the parity doubled states. This 
will show up in the present analysis in the meaning and quantitative estimate of 5% 
and the analogue to the Goldberger-Treiman relation. 

Thus, to better understand these issues it is interesting, if not essential, to study 
simple, solvable, strongly coupled toy field-theoretic models in which both heavy 
quark and chiral symmetries are present at the fundamental quark level, and the 
dynamics of chiral symmetry breaking is made explicit. We consider presently the 
simplest such scheme. We emphasize at the outset that this toy model is unrealistic 
and is intended only to convey the schematics of QCD chiral dynamics in heavy- 
light mesons (although we will brazenly attempt a fit to data). The simple model we 
consider is based upon a local gluonic current-current interaction Lagrangian: 

* j 

where the sums extend over all of the fundamental fermion flavors, both heavy 
and light, and we sum over the octet color index A. We view eq.(l) as essentially a 
QCD-inspired generalization of the Nambu-Jona-Lasinio model. For small g eq.(l) 
corresponds to the low-energy perturbative interaction generated by the exchange of 
a "massive gluon" of mass A/y/2. We propose to study this model using the technique 
of the large-A expansion, or equivalently, the fermion bubble approximation, with a 
cut-off at A. The model is exactly solvable in leading order. 

Solving the theory in the leading large- N approximation is equivalent to factoriz- 
ing eq.(l) into auxiliary fields describing the composite pions and heavy-light mesons, 
at the scale A and integrating out the heavy and light quarks to generate the effective 
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Lagrangian at a scale \x < A. The hadrons in our model, both light and heavy, appear 
as dynamically generated boundstates. In the light quark, meson sector we recover 
the chiral quark model of Manohar and Georgi [13] (with g\ — 1). In the heavy me- 
son sector we produce various boundstates of the heavy quark and the light quarks, 
and the full effective Lagrangian of these heavy meson boundstates coupled to light 
mesons is determined. The effective Lagrangian is manifestly heavy quark-spin and 
chirally symmetric. 

We will make certain further simplifying assumptions, keeping only terms in the 
renormalized effective Lagrangian that are ~ 0(1) or ~ 0(/xln(A//x)/A), while drop- 
ping subleading terms ~ 0(/i/A). We emphasize that we have in mind, presently, 
a hierarchy of scales, \x < A << M, where M is the heavy quark mass scale. The 
momentum-space loop integrals will extend from /i to the cut-off A. We view in the 
context of the model A to be a physical scale below which the theory is nonpertur- 
bative in g, but above which an effective softening of the point-like approximation 
due to the perturbative 1/q 2 gluon propagator takes place. While it is tempting to 
identify A with ~ Aqcd, we would hope that A ~ 1 GeV emerges from a fit to the 
physical quantities derived in the model. In fact, the simplest attempt at a fit to fs 
and f n yields A ~ 1.35 GeV, and most of the light sector observables are obtained 
within a factor of two. /i is an infrared cut-off which we would like to identify with 
the scale of light constituent quark masses. 

In the unbroken chiral symmetry phase the model produces the necessary de- 
generate parity doubling of the threshold spectrum of heavy mesons. In addition 
to the usual pseudo-scalar and vector HL mesons (the B and B* mesons which 
form a (CT, l - ) heavy quark symmetry multiplet), there is necessarily a scalar and 
pseudo-vector HL meson boundstate generated, which is a consequence of the chiral 
symmetry. This is identified with the sj l = | p-wave radially excited mesons (in 
the D-system this is distinct from the observed (D2) states which are sj' — | , (see 
Ming-Lu et.al. [14]). Unfortunately, these states have not been observed and will be 
fairly broad resonances, but they may ultimately be detectable. Technically, in the 
symmetric phase we must hold ji fixed at a nonzero value to protect from infrared 
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singularities. 

While the HQ symmetry maintains the degeneracy within the (0~, 1 _ ) and (0 + , 1 + ) 
multiplets, unbroken chiral symmetry implies the degeneracy of the two multiplets 
themselves. As we vary the model's coupling parameter to dynamically induce the 
chiral symmetry breaking, the theory develops a mass gap. This leads to a calculable 
mass splitting, elevating the (0 + , 1 + ) HQ multiplet and depressing the (0~, 1~). The 
mass gap between the groundstate mesons and the resonances is constant in the M — > 
oo limit and is given essentially by ~ gf n . This is the analogue of the Goldberger- 
Treiman relation of the theory, and is probably more general than our specific toy 
model result. Moreover, as a general result of the parity doubling, the axial vector 
coupling constant g& is not necessarily expected to be close to unity. In fact, g\ tends 
to be small based upon our fit, 0.32 (see Appendix B(ii); it occurs here as a term 
of order ln(A//i)/A, which is subleading to 1). This is a prediction which is thus far 
consistent with the upper limit in processes like D* — > D + 7T, though a measurement 
of the full D* width is still lacking to date. In the limit of very low q 2 pion emission 
we can decouple the heavier parity doubling states to return to the effective chiral 
Lagrangian for the (0~, 1~) groundstate mesons of ref.[6]. There remain in the low 
energy effective Lagrangian potentially important effects of the heavy resonances in 
chiral perturbation theory [12]. 

Thus, a key result we find is that the chiral mass gap, and hence an analogue 
Goldberger-Treiman relation, refers to the splitting between parity conjugate heavy 
meson multiplets in a heavy-light meson theory, i.e., heavy meson chiral theory is a 
parity doubled implementation of chiral symmetry. There are many other issues of the 
applicability of the chiral theory and its consequences which the present analysis will 
attempt to address. Though not entirely realistic, the model is completely solvable 
for various observables. /# is determined in terms of the short distance cut-off on 
the theory and the Isgur-Wise function is computed. The Isgur-Wise function result 
in the present model involves issues of going beyond the chiral logs, which arise also 
in matching composite mesons onto QCD. We will discuss this issue which is related 
to consideration of reparameterization invariance [14-16]. 
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II. Toy Model Field Theory with Chiral and Heavy Quark Symmetry 

(i) The Light Quark Chiral Dynamics 

Let us write the effective Lagrangian in the light quark sector, including the 
current-current form of the light fermion interaction Lagrangian of eq.(l): 

C = - m q )i, - y^^—^—il) (2) 

For concreteness we will take ip = {u, d), \ A are color matrices, and in the limit that 
the quark mass matrix m q — > 0, we have an exact chiral SU(2)l x SU(2)r invariant 
Lagrangian. This is a single gluon exchange potential, generated by a fake, mas- 
sive gluon of mass A/\/2. We treat the physics on scales q 2 < A 2 using eq.(2), in a 
fermion bubble approximation, imposing a UV loop momentum cut-off of A. Well 
above the scale A we would imagine the potential to soften into a 1/q 2 perturbative 
gluon exchange, hence A plays the role of a matching scale between strong infra-red 
physics and weak ultraviolet QCD. Finally, the "theory" consists of integrating out 
the fermions down to an infrared scale fi, keeping induced terms of order A 2 , A and 
ln(A//x) (we will discard terms that are finite, thus subleading, in the infinite A limit 
as a simplifying approximation). This generates an effective Lagrangian of composite 
particles. This is our essential approximation to the infra-red strong coupling behav- 
ior of QCD, or the "brown muck" of heavy-light physics. Overall, this is certainly 
a drastic approximation. Truncating on dim = 6 operators is, in a sense, a pure 
s-wave approximation to QCD, and cannot dynamically confine the quarks and dis- 
carding the subleading terms will limit the quantitative reliability of the model (the 
model could easily be improved). The physical value of A is determined in principle 
by fitting to the derived phenomenological parameters. The theory will contain the 
dynamical chiral symmetry breaking, and will determine a chiral Lagrangian of the 
heavy-light system. 

Upon Fierz rearrangement it is seen that the interaction Lagrangian of eq.(2) 
contains the Nambu-Jona-Lasinio model. The subsequent analysis is standard. We 
can factorize eq.(2) into a Yukawa theory with a static auxiliary field S = \{a+m a T a ) 
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on the scale /j, ~ A and then integrate out the fermions to determine the effective 
Lagrangian at scales /i < A. The field E is 2 x 2 complex at this stage, which implies 
parity doubling of the n and the parity partner of a, the i] is also present. This 
analysis is summarized in Appendix B. 

The light sector effective Lagrangian at scales /i < A can be identified with a 
linear cr-model: 



£ r describes the renormalized composite light mesons. We have written the renor- 
malized effective Lagrangian, so that g = g/\^Z^. Z 2 = (g 2 N/l6iT 2 ) ln(A 2 //i 2 ) is the 
finite, induced wave-function renormalization constant of the £ field. 

A (Tr£j£ r ) 2 term could be included in eq.(3), though it is subleading in N c , 
and for 577(2) x SU(2) with (a, n) real this is equivalent to the quartic term we have 
included. The theory can be tuned by choosing sufficiently large coupling g to develop 
a chiral symmetry breaking condensate, thus generating a constituent quark mass. 
The chiral symmetry breaking lifts of the isovector, + (Im(7r)) states The Re(ir) 0~ 
pion, becomes the Nambu-Goldstone mode. In QCD the residual U(l) symmetry is 
broken by anomalies and the effects of instantons. This generates additional terms 
such as an extra 't Hooft determinant, det £ + h.c. term, which elevates lm(cr) = rj. 
Any additional necessary Wess-Zumino terms should be incorporated as well. 

Since the light sector dynamics is not our principal concern in the subsequent 
analysis, we will henceforth assume that the fields (cr, n a ) comprising £ are real, so 
£ henceforth contains only the 0~ 7i a isotriplet and the real a isosinglet. Therefore, 
eq.(3) becomes a linear version of the chiral quark model ala Georgi, Manohar, and 
Holdom [14]. Nonetheless, we can dynamically put the model either in a symmetric 
phase, m 2 > 0, by choosing g 2 N/An 2 < 1, or in a chiral symmetry breaking phase 
m 2 < with g 2 N/47r 2 > 1. The critical bare coupling corresponds to m 2 = as 
fio — > 0. For further discussion of the light quark sector see Appendix B. 



ip(i$ - m q )^ - gip L T, r i/; R - gip R T,^ L 

+ Tr(^£j^£ r ) - m 2 a Tr(£];£ r ) + «Tr(m ? £ r + h.c.) 

+ATr(£t£ r £t£ r ) 



(3) 
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(ii) The Heavy-Light Quark Dynamics 

Now we focus on the dynamics of mesons containing one light and one heavy quark. 
The model produces one boundstate per channel in the fermion bubble approximation. 
We can conveniently solve the theory by factorizing the heavy-light (HL) interaction 
into auxiliary background interpolating fields with Yukawa couplings to heavy and 
light quark vertices. The original four-fermion interaction is recovered when the 
auxiliary fields are integrated out. Upon integrating out the quarks on scales A to 
//, the auxiliary fields acquire induced kinetic terms on the scale \x and thus become 
dynamical heavy-light mesons ( "i?-mesons" ) . In this way we derive the effective 
Lagrangian for the HL mesons coupled to the dynamical pions. 

The heavy-light fermion sector interaction Lagrangian, together with the HQ 
kinetic term, involves the HL cross-term of eq.(l) and can be written as: 

c HL = Q(i$ - M)Q -%Q^Y Q ^ r Y^ (4) 

Here we may generally take Q = (t,b,c.) to be a multiplet of Nh heavy quarks, 
and M the heavy quark mass matrix. We will presently consider, however, just a 
single heavy flavor in the following discussion, g should be viewed as the effective 
coupling at the scale A in both the light sector and the HL sector of our model. (In 
a more detailed discussion one might wish to distinguish the coupling constant in the 
heavy-light effective action from that of the light-light action at A; we ignore this 
possibility in the present paper). 

Upon Fierz-rearrangment of the interaction, again keeping only leading terms in 
1/Nq and writing in terms of color singlet densities, eq.(4) takes the form: 

Chl = Q( l ^M)Q + ^(WQ fl -^7W7U 

—jfTllditf'fQa ~ ^lulbtpi^lnlbQa^j (5) 

where 1 are the isospin indices, and a the heavy flavor indices. 

In the heavy quark limit we introduce a projection onto a heavy quark field with 
a well defined four-velocity v^. Presently we rewrite the full theory identically in 
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terms of a single four-velocity sector, corresponding to the four velocity of the heavy 
constituent quark or equivalently the boundstate heavy mesons: 

1 + 1) 

Q -> ~y~ exp(-iMv ■ x)Q{x) v (6) 

Note that (1 + f)Q v /2 = Q v , i.e., the field Q v always carries an implicit factor of 
(1 + f )/2. The HQ kinetic term then takes the form: 

Q>^Q m (7) 

The Isgur-Wise flavor symmetry is just the group of SU(Nh) rotations on Q®, and 
is now a manifest symmetry of our Lagrangian. We will consider just a single heavy 
flavor in the following. 

We now rewrite the terms of eq.(5) in a manifestly heavy spin symmetric form, 
letting Q — > Q v and further rearranging 7-matrices. Then, eq. (5) takes the form: 

g2 _. _. 

-QviJ-^^^-Y-i^Qv + Qvii^-^i^i^i^-Y-i^Q^j (8) 

We have now brought the interaction to a form which can be factorized by introducing 
heavy static auxiliary fields, (B,B'). To do so we must introduce four independent 
fields, B (B 5 ) are + (0~) scalars, while (B^) are \~ (1 + ) vectors. These form a 
minimal complete set of auxiliary fields needed to factorize eq.(8) in the HQ limit. 
Eq.(8) then becomes: 

£hl = QJv^df+Qv + gQ v i\)~B? v + igQ^^B^ 

-2A 2 (5l 5 ^ + W v B iv ) + 2A 2 (^ 5 5? + B i :B m ,) (9) 



Upon integrating out the B fields in eq.(9) we reproduce eq.(8). (Note that the B 
fields do not yet have canonical dimension of heavy meson fields; see Appendix A). 
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Eq.(9) is a heavy-spin symmetric form. We can assemble the auxiliary fields into 
complex 4 multiplets under 0(4) = SU(2) h x SU{2) h where SU{2) h {SU{2){) is the 
little group of rotations on Q v (ip and gluons) which preserves v^. One heavy spin 4 
multiplet consists of the + scalar together with the abnormal parity (1 + ) vector as 
(B : B 5tl ) (the four-velocity label „, and isospin 1 indices are understood): 

B' = (i-fB + ^B**) ff=(i^B + 1 ^)^-^j (10) 

The other 4 multiplet consists of the usual 0~ scalar and a 1~ vector (B 5 ,B^): 

B = (i^B' + ^Bn B={i^B* + 1 ,B»)( l -^-\ (11) 



Under heavy spin 0(4) = SU{2) h x SU{2)i rotations the (B,B 5>1 ) mix analogously 
to (B 5 ,B^). Note that v^B^ = always. We have introduced the caligraphic B and 
B' with the explicit projection factors. Falk has previously written similar effective 
"superfields" for excited mesons in model independent analyses; he includes an extra 
factor of 7 5 (relative to us) in his writing of effective fields for the (0 + , 1 + ) multiplet 
in a model independent approach [12]; for us the field B' has overall odd parity while 
B is even. 

The factorized heavy-light interaction Lagrangian then takes the compact form: 
C HL = Q v iv' i d lt Q v + gQ v (-iW i r i +&)il>i + h.c. 

+A 2 [Ty(BB) + Tr (#'#')] (12) 

Notice that the combination —i^B' 1 + B % is coupled. We emphasize that eq.(12) is 
exactly equivalent to the full four-fermion theory in the heavy quark symmetric and 
leading large- iV limit eq.(4). The theory forces a parity doubling of the heavy mesons 
upon us because the chiral symmmetry is controlled dynamically by g. For weak g 
the linear chiral invariance is realized and the theory must contain parity doubled 
meson states. Heavy spin symmetry organizes the parity partners into heavy spin 
4-multiplets. The effect of chiral symmetry breaking on the spectrum can now be 
investigated by solving the theory and choosing the broken phase. 
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See Appendix A for a discussion of normalization conventions. 



III. Full Effective Lagrangian 



We now proceed to "solve" the theory. The full effective Lagrangian for the 
heavy mesons is derived by integrating out the heavy and light quarks in eq.(12) 
over momentum scales fi < k < A, keeping the leading induced terms. Details 
of the explicit calculations are given in Appendix A. We begin the discussion with 
the use of the linearly realized chiral symmetry form, E = |(<7 + in ■ r), and we 
derive the nonlinear realization subsequently below. The loop integrations result in 
an unrenormalized effective Lagrangian. By performing a conventional wave-function 
renormalization and several field redefinitions we arrive at the full effective action valid 
to 0(/i/A) 2 : 

C LH = -i\ Tr(F; • dB) - i\ Tr(Fw • dB') 

Z z 



+ 
+ 



9r 

4 

h r 

4A 



Tr(Fx£>) - Ti(WaB') + Tr(Fvr • tB) + Tr(Fr • tB' 
Tr(F(a 2 + tt 2 )£) + Tr(F(5 2 + v 2 )B') 

Tr [B7 5 (0 71 ■ t)B - F7 5 {$ 7T • r)B' -W-f 5 ($a)B - B7 5 ($a)B' 



4/. 

+ A [Tr (BB) + Tr(Ftf') 



(13) 



The light quark PCAC masses are contained in the "shifted" a field, a = a + 
1m q \fZ^,l g. The parameters of this Lagrangian are determined as: 



2g 2 ^A 



V ^2 ^1 

A = i-(A 2 -Z 1 (A + /i)/2 7 r) 



z 1 



where: 



Z 1 = 



8vr 



(A-/i); 



9 2 N 
16tt 2 



ln(A 2 //i 2 



(14) 



(15) 
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The parameters defined above arise from the loop calculations presented in Appendix 
A. The g r , h r and k r are dimensionless. They are determined in principle by fitting 
the observables of the model as in Section IV. We will generally take /x to be of 
order the light quark constituent mass, and it will henceforth be neglected in the 
expression for Z\. Note that terms like B'^ 5 (cr$ a)B are potentially induced, but they 
are subleading as ~ 0(l/ln(A//x), relative to the terms we keep. 

We now identify the chiral representations of the composite fields in the effec- 
tive theory. This can easily be done by returning to eq.(12) and examining which 
heavy meson linear combinations couple to ipL and tpR. If we define the following 
combinations: 

1 .„„ „ 1 



B 1 



V2 



(B - iff) B 2 = ^(B + iff) 
v2 



(16) 



inspection of eq.(12) reveals that B\ (B 2 ) couples to ipR = (1 + / -f 5 )ip/2 (ipL — (1 — 
•y 5 )ip/2). Thus, the chiral representation of B\ (B 2 ) must be (0, \) ((^,0)). Writing 
in terms of £ = |(<r + in • r) (we will henceforth ignore the m q contribution which 
can easily be restored by shifting a — > a), the effective Lagrangian becomes: 



'LH — 



-i X - 1x{B x v ■ dB x ) - i X - Tr(B 2 v ■ dB 2 ) 



\Tr(BiT!B2) + Tr(B 2 T,Bi) 



h, 



2U 



Tr(BiBi) + Tr(B 2 B 2 ) 



(17) 



Inspection of the effective Lagrangian (as well as eq.(12)) confirms that it is manifestly 
invariant under SU (2) x SU(2) provided the fields transform as: 



(o,i) 



(18) 



We now see that indeed, eqs.(13, 17) have a structure analogous to that of a parity 
doubled nucleon theory, with B ~ (n, p)p =+ i, the normal even parity nucleon doublet, 
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and B' ~ {n,p)p=-i the odd parity doubling partner. We give a brief synopsis 
of such a system at the end of Appendix B. The essential results are that the axial 
vector current couples only through the perturbative k r term and describes transitions 
between parity partners, and the parity degeneracy will be lifted by (a). 

Note that eqs.(13, 17) describe the heavy meson dynamics in either a broken or 
an unbroken phase, i.e., it is simply a linear cr-model form. In the spontaneously 
broken phase of the heavy meson theory we can pass to the the nonlinear realization 
by replacing E with a unitary matrix field which is a function of angular pion fields, 
and a is now decoupled. Thus, the nonlinear realization is: 



£ = 2^ ex P( i7r ' T l '/*) 



2S// W = e (19) 
We can pass to the current form by performing the chiral field redefinitions: 

B 1 - ^B 1 £ 2 - ££ 2 (20) 

We then have the Lagrangian: 

9rf-K 



- l - Ti{B lV ■ {id + J L )Bi) - \ ^{B 2 v ■ {id + J R )B 2 ) 



Tr(fiifi 2 ) + Tr{B 2 B 1 ) + h.c. 
+A r [Tr(^i) + Tr(M 2 ) 

r Tr(^ 7 5 7^ 2 ) + Tr(B 2 ~7 5 7 ^i) 



where: 



(21) 
(22) 



and the chiral currents are: 



- -{J^r + J^l) 



(23) 

As usual the J a are matrices acting on the isospin indices of meson fields. The mass 
matrix of the chirally redefined heavy mesons is at this stage non-diagonal. We should 
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mention that if an extra 7 5 were included in the definition of the parity partner, then 
the axial current components of the J^l and J^r terms would carry 7 s factors, while 
no 7 5 would occur in the k r term. 

Note that the fields B\ and B2 are of mixed parity. The mass matrix can readily 
be diagonalized now that the Lagrangian is written in the current form: 

B = ^=(B 1 + B 2 ) & = -j=(B 1 -B 2 ) (24) 

with eigenvalues 2A r — gf n /2 and 2A r + gf n /2 respectively (recall that our normal- 
ization conventions imply the physical mass shift is 8M if the Lagrangian contains 
^5M(tvBB)] see Appendix A). The mass eigenfields are nontrivial functionals of the 
pions through the absorbed £, £' factors as in eq.(20). The Lagrangian now becomes: 

Clh = -~Tr(Bv-(id + V)B)-~Tr(Bv-(id + V)B') 
+ U r - 9 -^\ TtBB + (A r + ^\ TtBB' 

-- Tr(B(v ■ A)B') - - Tt{B{v ■ A)B) 

-tj TtB^^B + TrB^JjiB' (25) 

Note the appearance of the off-diagonal pionic transition terms of the form B(v -A)B'. 
[] At this stage it can be seen that these terms are associated with a Goldberger- 
Treiman relation, by taking = d^n ■ r/f n , integrating by parts, and using the 
equations of motion. One finds that the B'Btt amplitude has a coupling strength 
Qbb'-k = g r , and this is seen to be given by AM/ f n . 

We can decouple the heavier field B' to leading order in the mass gap g r f n by 
"integrating it out" (which amounts to setting it to zero in leading order). We can then 

1 The coefficient of this term corresponds to ft = 1 and k r — g in Falk's notation [12]. Our 
conclusion is that k r = gA < h = 1, and following Falk's analysis the chiral perturbative contribution 
of these resonances, e.g., to fDs/fo, is significant. 
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perform the residual mass redefinition: B — > exp(— iMv -x)B where M = 2A r —g r f 7T /2 
to yield the final result: 

C LH = -\ Tr Biv ■ (d + V)B-i^-TrB-f 5 4 B (26) 
where we now discover that: 

g A = k r = 2gf n y/z~ 2 /Z 1 (27) 

Our fit to the model yields g A = 0.32 (see eq.(38) and discussion). Eq.(27) is equiva- 
lent to the point of departure taken by ref. [6] in writing effective Lagrangians involving 
simultaneous heavy symmetries and chiral symmetries. Use of this effective Lagragian 
is justified so long as q is small compared to the mass gap. We see that g A here arises 
from the perturbative k r term, which is subleading to unity in our expansion. 

In summary, the central observation of this analysis is that the underlying chiral 
representations of the full HL meson theory is a parity doubled scheme. There are two 
general implications of such a scheme: (1) The mass gap between the parity partners 
arises from VEVa. Thus a Goldberger-Treiman relation refers, not to the overall 
mass of the B mesons ~ M, but rather to the mass splitting between the even and 
odd parity multiplets: 

AM = gBB'irfir (28) 

Here gBB'n = flV is the BB'n transition coupling constant and is the analogue of the 
gNN-w in the nucleon system. We note that the light quark constituent mass is given 
by m c ~ grfn/2 so we expect AM 600 MeV, however this must be obtained in 
principle from a fit of the model to all data (see section IV. (ii); Unfortunately, without 
exceptional circumstances the width of this state is too large for observation.) (2) 
g A is not necessarily expected to be ~ 1, being given by a subleading perturbative 
contribution, k r , alone. This is essentially a consequence of parity doubling and 
contrasts the or the chiral quark model in which, g\ — 1 is a leading term. The fit 
we present below in section IV(ii), which is crude, yields g A 0.32. This result may 
be indicated in the D-system where D* — > D + it gives g A < 0.7 [6, 11]. 

IV. Other Observables: f B , and Isgur-Wise function 
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(i) Heavy Meson Decay Constant, f B 



We presently compute the heavy meson decay constant f B . Consider the heavy- 
light axial current Qj^ip. We can compute the renormalized matrix element: 



-i 



Z 1 J d'x e-* M ™ (B\ Qixh^ix) |0) = f B ^M B v, + ... 



-iMv-x 



(29) 



As a consequence of the heavy quark spin symmetry, £> 5 and have identical decay 
constants for the axial vector current, while B and B^ have the same decay con- 
stants for the vector current. The _B-meson must have a properly normalized kinetic 
term, which includes the finite renormalization effects, B — > \J Z{ l B. We adopt a 
conventional normalization in which we expect f B ~ 180 MeV. 



The amplitude on the Ihs takes the form: 
igN f d 4 k 



I (2^ ^IVA* - 1 + m q )(l ~ f ) (^7 5 5 5 + 1 V B V )\ 



1 



2^Z\J (2tt) 4 



((A;-p) 2 -fi 2 ) (oi) 



1 



A 2 — /i 2 + 7Ti> • p(A - /x) + -ttcFA - f2 2 ln(A 2 //i 2 ) + • p) 2 ) 



(30) 



We see that the integral involved here is identical to 1\ of eqs.(60, 62), and thus the 
equations of motion can be used for the B 5 fields. Upon use of the equation of motion, 
shifting v ■ p — > 2 A + ... = 2A 2 /Z 1 — A/ir + ... a large cancellation is seen to occur on 
the rhs of eq.(30) leaving: 



2gN 



2%A 3 /Z 1 



A'/JZ 1 



We thus obtain: 



{MbJb = (A) 3/2 



4V27T 



(31) 



(32) 



gWN 

For example, let g 2 N/4ir 2 = 1 and use f B = 180 MeV, M B = 5 GeV as input 
parameters, to find A = 1.35 GeV for Nc = 3. Remarkably, our result is insensitive 
to the light quark masses. 
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fs is a measure of the wave-function of the meson at the origin in a nonrelativistic 
potential model. We can compute the wave-function in principle in our model by 
point-splitting the current in eq.(29): 

y/zT* J d 4 x e~ lMv - x (B\ Q(x - e/2)i) ^{x + e/2) |0) = A/M>(e) (33) 

where J\f is a normalization factor. This essentially replaces the momentum space cut- 
off procedure by a point-split regulator, and A ~ 1/e. The wave-function is singular 
at the origin, and is not normalizable without a spatial cut-off of the normalization 
integral at 1/A (our theory makes no sense at shorter distances than this). Thus the 
wave-function is at the origin is given effectively by: 

|^(0)| ~ K 2 / \Jz\ ~ (A) 3/2 . (34) 

This implies that the result for /# is insensitive to infrared parameters such as the 
light quark masses in our model, and indeed we find Jb, s — J'b,u- This is a defect of the 
model, but it is an expected result of an extremely relativistic, potential dominated 
system. In this sense, QCD lies somewhere between this extreme result and that of 
a naive potential model. 

(ii) Fitting the model to data: 



While the model we have presented is not likely to be quantitatively successful, 
we can attempt a fit to observables, and predict some features of the HL meson 
system. We use as independent inputs f w = 95 MeV, and fg — 180 MeV for M B = 5 
GeV. The latter implies A = 1.35 GeV as discussed in the previous subsection. We 
see, owing to the smallness of the ratio f 2 /A 2 = (k — l)/g 2 , that k = g 2 N/An 2 is 
very close to unity. In defining Z 2 we cut-off the renormalization group flow at an 
infrared scale \i ~ m c taken as the approximate constituent light quark mass. Then, 
to obtain Z 2 = /m 2 c we self-consistently solve for the constituent quark mass 

m c = \gfn/ ^Z 2 {m c ). This yields: 

q 2 N 

^-=- = 1.065; # = 3.75 A = 1.35 GeV Z 2 = 1.11 m c = 169 MeV 

An 2 , 

(35) 
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m c is about a factor of two too small. We can moreover use the pion mass, m n , to 
extract the light quark PCAC masses: 

mlfl = cu(m u + m d )U ^ = /f^- = 0.25 {GeVf (36) 



87T 2 v^ 2 

hence, m u + = 8.6 MeV, which is to be compared with the conventional ~ 15 
MeV, and is small. Also, m s w m 2 K (m u + rridj/m 2 ~ 107 MeV is small. 

The mass gap between the excited + and groundstate 0" mesons is then: 

AM = g r f n « 2m c w 338 MeV [600 MeV] (37) 

The result in brackets obtains when the known constituent masses are inputted. The 
decay width T(0 + — > 0~1) is given by (AM/ f n ) 2 \k n \/8iT. This is much too large for 
observation of these resonances when ~ AM ~ 600 MeV; with the lower estimate 
of AM ~ 338 MeV the width approaches 150 MeV, which is still too large. Hence, 
the direct observation of the parity partners of the groundstate is unlikely. Their 
effect in chiral perturbation theory is nontrivial [12]; conceiveably the decay width 
T(D S (1 + ) — > D^ d (l~) + K) ~ \p K \ is kinematically suppressed by the if-meson mass 
and the 1/M corrections to the D masses. 

We obtain the axial coupling constant: 

2^^ = 4^1np^^ 032 



We might expect both AM and g& to be underestimated in this approximation, as 
are the light sector observables, owing largely to the short-distance singularity of our 
wave-function. 

qa can be in principle extracted from the decay D* + — > D° + 7r + , though it is un- 
measured to date. This decay partial width is given by ref . [6] , and in our conventions 
it takes the form: 

91 VI 3 (39) 



where ~ 38.9 Mev. While this width is not yet measured directly, we can use the 
analysis of ref. [19] to obtain an estimated result of T = 53.4 KeV from the measured 
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branching ratio of (D* — > D^)/(D* — > L>7r) and a potential model calculation of 
Z)* — > D7. Combining these results we find g\ = 0.56, which is compatible with 
the parity doubled interpretation, but is also not far from the naive g& ~ g\ ~ 0.8 
from the chiral quark model (note that we derive the chiral light-quark model here 
with g\ — 1, thus our prediction of g\ ~ 0.3 represents a significant suppression). 
Amundson et.al. [11] give the current experimental limit of gA < 0.7 consistent with 
this result. Thus, our model indicates that g^ is suppressed and smaller than unity, 
giving the physical underlying rationale, though the situation is arguably not decisive. 

Note that Z x = ttA/2 ~ 2.12 GeV and Z 2 ~ 1.1. Hence, 2 A = 3A/?r « 1.3 GeV. 
Our model seems to suffer from generating a value of A that is slightly large. This 
implies e = 4Z 2 A/Zi ps 1.2, suggesting that our approximation of truncating on the 
Z 2 (v ■ p) 2 /Zi terms is probably unreliable (Appendix A). 

The binding energy is determined in the model. Neglecting the light quark PCAC 
masses we have in the infinite mass HQ limit: 



where m c is the constituent light quark mass (the latter term is small, but non- 
negligible). For the fit we have presented we find h r = 2g 2 ^fZ~ 2 kjZx = 17.9, and 
g r = gj\[Z~2 = 3.55. If we use a conventional charm (b-) quark mass of 1.2 — 1.8 
GeV, (4.5-5.0 GeV) this overestimates: M D ^ 2.4-3.0 GeV, ( M B = 5.8-6.3 GeV). 
These results should be corrected for finite mass of the heavy quark. The corrected 
boundstate mass is: 



This yields a result M D w 2.0 - 2.75 GeV, ( M B = 5.65 - 6.3 GeV). This illustrates 
the problem of A being too large in the model. 

The effect of the explicit SU(3) breaking light quark masses is calculable, upon 
restoring these terms in eq.(13) as contained in the shifted a field. Using the full 
constituent quark mass m c = g r a/2 we have: 



M D)B = M cf) + 5m; 



5m = 2 A — m c + 




(40) 




(41) 




20 
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(m c - m d ) + (2.1 x 10~ 3 {MeV)- l ){m 2 c - m 2 c ,) 



(42) 



For the D s -Dq mass difference we take m c = 450 MeV (the strange quark con- 
stituent mass) and m c > = 300 MeV to obtain Mo s — Md(u,<i) — 86.25 MeV. (If we 
use the predicted m c r = 169 MeV and m c = 276 MeV we obtain 52.5 MeV). This 
compares to ~ 100 MeV experimentally. It shows, however, that the model must 
include the effects of the a 2 term in computing these differences. The Mo+ — M^o — 
[(+0.26), (— 0.3)] (m^ — m u ) ~ (2.6, —3) MeV (using standard constituent masses in 
the first entires, and the model's derived constituent masses in the second). This is 
subject to electromagnetic corrections, estimated to be +2.0 MeV. 

We have seen that fs is insensitive to the light quark masses in this model. Thus, 
we obtain f Bs / fsu^d = 1) while lattice results yield ~ 1.09 [20]. This result owes to the 
unrealistic non-normalizeable singularity of the wave-function at the origin. This is 
consistent with the behavior of the binding energy for small constituent quark mass, 
in which increasing the constituent mass actually decreases the meson mass (for large 
constituent mass the a 2 terms contribute to increase the meson mass). 

(in) Isgur-Wise Function 

The analysis of the Isgur-Wise function in the model involves a careful treatment 
of the cut-off procedure. We select a preferred cut-off by demanding the validity of 
reparametrization invariance (or the residual mass symmetry) [15-17]. 

We consider the transition amplitude in 4-velocity defined by the matrix element 
(B v \ Q V TQ V > \B v i), where T is an arbitrary Dirac matrix. 




(43) 



This involves the integral: 



h 




A-(v^ + v' 11 ) 



(44) 
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where the latter term follows by symmetry, since A can only depend upon v ■ v' and 
is even under v <-> v'. Now multiply by v + v'\ 



2A ■ (1 + v ■ v') 



Therefore: 



d 4 k k-(v + v') 1 1 
(2tt) 4 (k 2 -m 2 ) vkv'-k 

d A k 2 1 
(2tt) 4 (k 2 -m 2 )vk 



— (A - 2m) (45) 



(A - 2m) (46) 



167r(l + v ■ v' 

and we conclude that the Isgur-Wise function is given by: 

£(v ■ v') = (47) 

This should be true if the momentum space integral is Lorentz-invariant and finite. 
Computing the integral directly, without recourse to the symmetry argument one can 
obtain: 

1 , A / >7r /2 cos# ,^ oX 

2^ Wo d V2W cos sin ^ (48) 
which agrees with the previous result. 

This result contains a t-channel pole at t — —(Mi + M 2 ) 2 , where M 1 (M 2 ) is the 
incoming (outgoing) heavy meson mass. One might ask if this is consistent with the 
slope constraint of de Rafael and Taron [18] arising from t-channel unitarity? Our 
slope, £'(0) = —1 is inconsistent with the their lower bound of —1/2 arising from a 
t-channel branch-cut at threshold. Grinstein and Mende [18] have pointed out that 
the de Raphael-Taron constraint is weakened by effects of resonance poles, as we 
are presently observing. However, the t-channel unitarity constraint is an interesting 
issue in HQET. In an HQET such as we have studied, the anti-particle has been 
discarded at the outset, and with it goes crossing symmetry and t-channel unitarity. 
Moreover, our cut-off theory would seem to require the bound of Q 2 < A 2 without a 
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unitarization. Since Q 2 = 2M 2 (1 — v- v'), we see that this bound corresponds to the 
limit v ■ v' — > 1 for M — > oo. Nevertheless we can compute the t-channel behavior 
by incorporating the heavy anti-quarks and computing the large- N bubble sum with 
the full interaction. While we do not present this analysis here, we find, perhaps not 
surprisingly, a Nambu-Jona-Lasinio pole at Mi + M 2 is generated, and our slope is 
consistent with the existence of this pole. 

The previous result of eq.(47) is, however, sensitive to the definition of the cut-off 
procedure, which we have taken to be a Lorentz-invariant Euclidean momentum space 
cut-off. Different results follow if the energy integrals are first performed by residues, 
and then a 3-momentum cut-off procedure is used. To see this let us compute directly 
with a 3-momentum cut-off Let v' = (1,0) and v = (vq,v). First we perform the 
energy integral by closing dko below to pick up the single pole: 

f d 4 k I k + v k - 

A 1 + V • V ) = / — =; 

1 ' J (2vr) 4 ^ kl - k 2 + 

= A / dk 

4tt 2 Jo 

(49) 

This result, using a non-Lorentz invariant regularization procedure, differs signifi- 
cantly from eq. (47) in which the Lorentz invariant cut-off was used. 

There is however an implicit gauge invariance in heavy quark effective theories 
associated with the "residual mass ambiguity." One is free to add a term xQ v Qv to the 
effective Lagrangian of eq.(12). x should be viewed as a gauge potential in the sense 
that if we redefine the heavy quark mass M — > M+/i, and thus Q — > exp(i/iv-x)Q we 
can compensate this gauge transformation by shifting x —> X + 1 1 - Hence iv ■ d + x is a 
covariant derivative. This is essentially the demand that the global zero of energy of 
a classical theory be arbitrary. This symmetry and its implications will be discussed 
elsewhere, however we can see immediate implications for our present problem. 

We can observe that the non-Lorentz invariant regularization procedure violates 



1 



,i>o&;o — k ■ v 




In 



v ■ v') 2 — 1 



v ■ v' + y (v ■ v') 2 — 1 
v ■ v' — \ (v • v') 2 — 1 



+ 2 
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the x symmetry. Consider the integral involved in our calculation of the Isgur-Wise 
function: 

f d 4 k 1 
~ J (2tt) 4 \{k + p) 2 - m 2 }{v k + n) (5 °^ 
We have chosen to route the external momentum p through the light fermion line. 
The x symmetry applies to the external heavy mesons and requires that the following 
shift in v ■ p and fi be a symmetry of the integral: 

vp^v-p + x; /i^^ + x (51) 

This is readily seen to be a symmetry in the case of the momentum p routed through 
the heavy fermion line. 

In the present example we can implement this by shifting p — > p + v x- Therefore 
the shift in the integral is: 

5 r d 4 k { -2vk 1 \ 

5x~ ~ J (2tt) 4 \ [{k + p) 2 -m 2 ] 2 {vk + fi) ~ [{k + p) 2 - m 2 ]{v k + /i) 2 J ( ' 

The symmetry condition is -^1 = 0, and is equivalent to demanding that the integral 
generates no nontrivial surface term upon shifting k — > k + a. For simplicity we 
consider the surface term: 

r ^k_ ( —2v ■ k 1 \ , , 

J (2vr) 4 \{k 2 -m 2 ) 2 {vk) (k 2 - m 2 )(v k) 2 J { ' 

If we evaluate S using the covariant cut-off we find that 5 = 0. 

Now consider computing S by first performing dko by residues, then the residual 
3-momentum integration with a cut-off. We find: 

d 4 k v ■ k f d 4 k 1 



(2tt) 4 (k 2 - m 2 ) 2 (v ■ k) J (2tt) 4 (k 2 - m 2 ) 2 

k 2 dk 



47r 2 Jo (k 2 + m 2 ) 3 / 2 



(54) 



Consider now 



d 4 k 1 —i r in r°° k 2 dk 

sin OdU 



Jo 1 



(27r) 4 (A; 2 -m 2 )(t;-A;) 2 8tt 2 Jo Jo (k 2 + m 2 ) 3 / 2 (v - \v\u(k) cos 9) 

(55) 
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where u(k) = k/^/k 2 + m 2 and = (v ,v), and thus v% — v 2 = 1. Notice that if 
either m — > or if {7 — > then 5 — > 0. Let us expand in m 2 , using the latter results, 
to find for S: 

_j_ _ 2 f°° x 2 dx 
b ~ An 2[i Vo) Jo (l + x 2 )^ 

= I^ 1 "^) 2 ) ( 56 ) 

Here we introduce a four- vector = (1,0) which is the direction of the dko line 
integration. 

This latter result implies that the \ symmetry is broken when the ko line integral is 
not parallel to v^. For the computation of the Isgur-Wise function where v ^ v' then 
the x symmetry can never be present in the residue computation. However, utilizing 
the Lorentz invariant cut-off we see that the x symmetry can be maintained. The 
X symmetry therefore requires that we reject the result of eq.(49) in favor of eq.(47) 
which is consistent with the absence of momentum space surface terms, and the 
attendant symmetry. 

VI. Conclusions 

We have presented perhaps the simplest, solvable, strongly coupled toy field- 
theoretic model in which both heavy quark and chiral symmetries are present at 
the fundamental quark level, and the dynamics of chiral symmetry breaking is made 
explicit. We find that the chiral representations of the heavy mesons are parity 
doubled. This has a well defined meaning in the toy model because we can tune 
the coupling constant to restore the spontaneously broken chiral symmetry. In the 
symmetry limit the groundstate is a degenerate system of (0~, 1~) and (0 + , 1 + ) heavy 
mesons. When chiral symmetry is broken the degeneracy is lifted, elevating the 
(0 + , 1 + ) and depressing the (0 - , 1~) heavy meson multiplets. We obtain the full 
chiral Lagrangian containing the parity doubled composite HL mesons together with 
the composite pions. The mass gap between the multiplets is given by gf n , and 
the analogue of the Goldberger-Treiman relation of the system reflects this, Qbb'-k = 
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AM/f w . We are able in the broken phase to pass to a nonlinearly realized chiral 
symmetry, and to write a purely derivatively coupled pion effective Lagrangian. We 
can then decouple the heavier parity doubling field to arrive at the conventional low 
energy effective chiral Lagrangian for B. 

We believe that the general phenomenon of the parity doubled chiral representa- 
tions of heavy mesons is inherent to QCD itself. We emphasize at the outset that this 
toy model is intended to convey the schematics of QCD chiral dynamics in heavy-light 
mesons The model is designed to imitate these dynamical features of QCD, rather 
than provide a detailed phenomenological fit. Nonetheless, the simplest fit seems 
to agree within a factor of two to the expected values of physical quantities, and is 
predictive. While we would be inclined to trust the result qa = 0.32 only to within a 
factor of 2, the model suggests that qa is smaller than might be naively expected on 
the basis of the simple constituent quark model in which. The direct observation of 
the parity partners of the groundstate mesons is unlikely owing to their large widths. 
It would be interesting to extend these results to the heavy quark containing baryons 
where similar conclusions must hold. 

Our analysis achieves the basic systematics of chiral symmetry in these systems 
where we might expect potential models to fail. The chiral symmetry limit is relativis- 
ts, and the chiral symmetry breaking is a dynamical rearrangement of the vacuum, 
two features which would be hard to realize in any potential model treatment. One 
must be careful in estimating the value of qa in a naive potential model unless the mix- 
ing with the parity doubled states is under control. As we have observed in eq.(21), 
the qa term is a transition matrix element between the + and CT states in the mixed 
parity basis appearing there. In a basis in which the qa term is diagonal, the mass 
matrix must be correspondingly diagonal. There remains the transition amplitude 
term between the parity partners (some authors include an extra factor of 7 5 in the 
odd parity fields, and this transition term can then be mistaken for the qa term in a 
mixed parity basis). In our model, the constituent quarks are found to have g\ — 1, 
and yet the value of qa obtained in the Lagrangian of ref.[6] is suppressed to ~ 0.32. 
This is a subtlety of parity doubling which must be treated with some care. The res- 
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onances may have important contributions in chiral perturbation theory to quantities 
such as fosl Id and flavor ratios of Isgur-Wise functions [12] (in the notation of Falk, 
h = 1 and qa = g, and thus h 2 » g 2 in our model, so the resonance contributions 
are significant). 

We have studied the physical predictions of this system. The wave-function of 
the theory is too singular at the origin to represent a realistic QCD wave-function. 
This is a consequence of the strong coupling of the point-like four-fermion interaction 
term. While it is a defect of the model, it indicates the trend in a theory in which the 
potential term is dominating the dynamics. For example, we obtain the unrealistic 
fBs/fBu,d = 1, because the singular short-distance behavior of the wave-function 
becomes insensitive to the infra-red parameters of the theory. This contrasts lattice 
results, indicating fss/ ' fBu,d = 1-1 [20]. However, a weakly coupled potential model 



In our analysis we fix A ~ 1.35 GeV from f B and examine the relationship with 
the cut-off wave-function at the origin. Inputting also /„- fixes g, and marginal results 
(within a factor of 2) obtain for AM, g A and the light quark sector. A defect, related 
to the short-distance singularity of the wave-function, is the fact that for small light 
quark constituent mass, the groundstate mass is actually depressed as the light quark 
constituent mass is increased from small constituent mass. Nonetheless, the common 
h r term is sufficiently large for m c ~ 300 MeV that a reasonable result for M B s — M B u,d 
emerges from the fit. 

Of further interest is the Isgur-Wise function, which is associated with an am- 
biguous linearly divergent integral in the present scheme. The ambiguity is resolved 
by invoking "residual mass invariance" [15,16], or equivalently, "reparameterization 
invariance," and enforcing an associated Ward identity [17]. The simple Isgur-Wise 
function corresponds to a t-channel threshold pole at (Mi + M 2 ) 2 . This pole is beyond 
the cut-off scale of our model, but it may be indicative of a Nambu-Jona-Lasinio re- 
sult when the QQ system is studied. In fact, the fundamental issues raised by de 
Rafael and Taron can in principle be explored in this scheme [18]. We will defer this 
discussion to another place. 



would give the larger result Jm s /md ~ 1. 




2 [21]. 
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We believe there remains much to do in dynamical analyses of this kind for heavy- 
light systems and their interactions. Our model has inherent shortcomings. While the 
agreement of this crude model with observation is marginal at best, it suggests that 
improvements, such as a Pagels-Stokar approximation, Holdom's approach [22, 23], 
or "Russian sum-rule" methods [24], will lead to more reliable estimates of crucial 
heavy meson observables. The singular behavior of the wave-function is not expected 
in a more realistic scheme. Replacing our pure s-wave dynamics by QCD ladder 
approximation is clearly of some interest. For example, pinning down a prediction 
of qa or the Isgur-Wise function from such models would be quite interesting. The 
full range of phenomenological applications of generalized models would seem to be 
an interesting direction for future research. This toy scheme is a first step in that 
direction and highlights the challenges and advantages for more elaborate approaches. 
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Appendix A: Fermion Loop Approximation 



(i) Zero momentum pions 



Let us now integrate out the heavy and light fermion fields in eq.(12) to pro- 
duce an effective Lagrangian for B and B' . This can be viewed as a block-spin 
renormalization of the theory of eq.(12) defined at the scale ji — A, to a new scale 
fj, < A, and is analogous to the treatment of the light quark dynamics in Appendix 
B. We begin in the approximation of treating the a and n fields as zero-momentum 
(constant in spacetime) backgrounds (small momentum n amplitudes are considered 
subsequently). We note that the fermion propagators take the form: 



S HQ {k) 



v ■ k 



S LQ {k) 



. ( k 1 +m q + gZ 5 
1 \ k 2 -9? 



(57) 



where: 



Q- = (,„„ + (jo I2f + <j 2 n 2 / I: K 2 = = -a + ■ r (58) 



We obtain from the diagram of Fig.(l) (recall that the B contain (1 +^ )/2 projection 
factors) : 



iS BB = ~9 2 N J 



d A k 



Tr 



(k- P ) 2 -n 2 / 

= -g 2 N Tr [(-iF 7 5 + B)h{-i^B' + B) 



v ■ k 
(59) 



and: 



d 4 k ({k 1 -i))+m q + gZ 5 



(k — p) 2 — il 2 J v ■ k 



(60) 



We carry out a "block-spin" integration over heavy and light quark modes between 
the scales fi and A in Euclidean momentum space. The integrals are evaluated with 
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a Euclidean 4-momentum cut-off: 
d 4 k 1 



I 



(27r) 4 ((k - p) 2 - £l 2 )v ■ k IQtt 
d 4 k 1 i 



(2tt) 4 ((A; - p) 2 - il 2 ) 



IQtt 2 

V 



—l- 



(2vr) 4 (k 2 - n 2 )v ■ k IQtt 2 



A 2 - /i 2 - f] 2 ln(A 2 //i 2 ) 

A 2_ /U 2_ fi 2 ln(A 2 //i 2 ) 



d 4 k 



k,. 



ln(A 2 //. 2 ) 



(2tt) 4 (k 2 - il 2 ) 2 v ■ k IQtt 2 
Then Tr(7i) can be written as (note Tr(B-f 5 B) = Tt(B'B) = 0, etc.): 
ig 2 NTr(BI 1 B) 

= - l - Tr(BB)[(v ■ p + ga/2){Z x + AZ 2 v ■p) + Z 1 (A + /i)/tt - 2\Q\ 2 Z 2 ] 
ig 2 NTi(B'{-i^)h{-i^)B') 

= -\ Tt(B'B')[(v ■ p - ga/2)(Z! + AZ 2 v ■p) + Z 1 (A + /i)/vr - 2\tt\ 2 Z 2 ] 



(61) 



ig 2 N[Tr(B (-rf)hB) + Ti {Bhi-t^B')} 

= ~ Tt[B\ 9 tt ■ t/2)B + B(grT ■ r/2)B']{Z l + AZ 2 v ■ p) 



where we let go — go + 2m q and 

g 2 N 



Zx = 



Z 2 = 



S 2 iV 



[ln(A 2 //i 2 ) 



(62) 



(63) 



8tt v ~~ r ~ n * \Qtt 2 

(note that the expression for Z\ contains a factor of 1/tt, not 1/tt 2 ). 

(ii) The gA term 

Now consider small, but nonzero (a, tt) momentum q^, with q 2 pa 0. We compute 
the effective Lagrangian, where the (a, tt) are coupled through S 5 . We then have the 
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amplitude of Fig. (2): 



iS 



BBY, 



d A k 



Tr 



(2^) 4 
HF7 5 + B) 



(ff +^/2)(q + z7r-T 7 5 )(^ -rf/2) 
((A; + g/2) 2 - fi 2 )((A; - g/2) 2 - fi> ■ (k - p) 



-i^B' + £) 



(64) 



We are interested in the divergent terms of order q, since the q = term has previously 
been computed: 



— 9 3 N 



W) 4 



Tr 



[ft J](a - in ■ r-f 5 
(k) 2 (k) 2 (vk) 



i'fB' + B) 



~4^ 3jV 16^ Tr \^- lB '^ + B)[t,6K<r-™- ri 5 )(-ti 5 B' + B)} ln(A 2 //i 2 ) 



= --<?Z 2 Tr [H#' 7 5 + , jtf ](<r - in ■ t 7 5 )H 7 5 B' + B) 



(65) 



If we now expand the result of eq.(65) we observe some simplifications, e.g. B[f ,$]B 
0, and we obtain: 



-gZ 2 Tr [Bi {in ■ r 7 5 )B - B'fr {in ■ t^)B' + iB'^fr oB + Bi^i oB' 



(66) 



This implies an operator in the effective Lagrangian of the form: 

= X -gZ 2 (Tr [B^^B - W^^B'} d»n a - Tr [F 7 5 7 ^ + B 7 5 7 M #] d»o 



(67) 



(Hi) Normalization Conventions 
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Consider a complex scalar field $ with the Lagrangian: 

dp&dPQ — (M + 5M) 2 &$ (68) 
Define $' = \/2M exp(iMv ■ x)<& and the Lagrangian becomes to order 1/M: 

iv^d^' - 5M$' f $' (69) 
Now let $ = |(1 + f )«7 5 $' and write in terms of traces: 

- i\ Tr(&v ■ d$) + 5Afi Tr($ t $) (70) 

Thus when the Lagrangian is written in terms of B and B' the normal sign conventions 
are those of the vector mesons, and opposite those of scalars, i.e., the term in the 
Lagrangian +|<5MTr(£>£>) an increase in the B 5 mass by an amount 5M. A properly- 
normalized kinetic term is —i \ Tr(Bv ■ dB), with the overall minus sign and \. 

(iv) Structure of the Effective Lagrangian 

The heavy meson effective Lagrangian therefore takes the form: 



-i]-Z 1 Ti(Bv ■ dB) - i\z x Hx{B'v ■ dB') + 2Z 2 Hx{B(v • d) 2 B) + 2Z 2 Tr(B'(v ■ d) 2 B') 



g 2 z 2 



Tr{BaB) - Tr(B'aB')] - igZ 2 [Ti(Bav ■ dB) - Tr{B'av ■ dB') 
Tr(Fvr • tB) + Tr(B~7T ■ tB')] - igZ 2 [Tr(F7r • rv ■ dB) + Tr(B7r • tv ■ dB') 



Ti{B{d 2 + 7i 2 )B) + Tr{B'{a 2 + ti 2 )B') 



+ (A 2 - Zi(A + m)/2tt) [Tr(BB) + Ti(B'B" 
+^-gTr [B 7 5 ^ (tt • t)B - W~f 5 $ (tt • r)B' - F 7 5 ^ (tr)fi - B 7 5 ^ (tr)fi' 



(71) 
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If we define: 



T = [1 + {AZ 2 /Z l )ivd]- 112 



(72) 



then eq.(71) becomes more compactly: 



1 



C LH = -i-Z 1 Tr(TBv-dTB)-i-Z 1 Tr(TB'vdTB') 

Z z 



4 

gZi 

4 

9 2 Z 2 



Tr{TBaTB) - Tr {TB'aTB') 



Tr(TB'ir ■ tTB) + Tr(T#7r • tTB 1 
Tr(B(a 2 + n 2 )B) + Tr(F(cF 2 + n 2 )B') 



+ 



(A 2 - Zx(A + m)/2tt) [Tr(fifi) + Tr(F-B']) 



+-^0 Tr [i3 7 5 ^ (tt ■ r)B - i3' 7 5 ^ (tt • r)B' - B'^$ (a)B - &y 5 ? (a)B' 
z 

(73) 

To simplify the subsequent analysis we will assume that the subleading terms of order 
Z 2 v -p/Zi are negligible, and take T = 1. Since these terms arise upon expanding 
the loop integrals in powers of 1/A, we cannot self-consistently use the effective 
Lagrangian in this form unless this condition is at least approximately valid. We see 
that other terms, such as the last one in eq.(73) which leads to qa, are leading in this 
order and describe various physical processes. Thus, we expect the amplitudes these 
terms describe to be small. If 4Z 2 v ■ p/Zi is large, then we must retain full analytic 
expressions for the loop integrals to fit the theory. 

We see that there is thus an induced kinetic term for the B and B' fields with 

a common wave-function normalization. We absorb the factor Z\ into the fields 

-1. 



as B — > \[Z\ B. Thus, with the field redefinition we then have the full effective 
Lagrangian: 



C 



LH 



-i\ Tv{Bv ■ dB) 



i- Tr(B'v ■ dB 1 ) 



Tr (BaB) - Tr(B'aB') 
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Tr(B 7 7r • tB) + Tr(Bvr ■ tB') 



+ 



9 2 Z 2 



Tr(B(a 2 + tt 2 )B) + Tr(B'(a 2 + n 2 )B') 



AZ 1 

+ A [Ty(BB) + Tr(FB') 



2Zi 



g Tr i3 7 5 ^ (tt • r)B - B'f$ (vr • r)B' - B'f$ (a)B - Bf$ (a)B' 



where: 



A = ^(A 2 -Z 1 (A + /U )/2tt) 



(74) 



(75) 



The equation of motion in momentum space is v • p — 2 A + ... and 2 A is the mass 
difference between the heavy meson and the heavy quark in the chiral symmetric 
phase: 

M B = 2A + M Q (76) 
Note that A > (A < 0) for g 2 N/16n 2 < 1 (g 2 N/16n 2 > 1). 
Appendix B: Light Quark Dynamics 

(i) Deriving the Constituent Quark Model 



The effective Lagrangian in the light quark sector is: 

_ g 2 _ \ A _ \ A 



(77) 



For concreteness we will take ip = (u, d), and in the limit that the quark mass matrix 
m q — > 0, we have an exact chiral SU(2) x SU(2) invariant Lagrangian. This can be 
viewed as a single gluon exchange potential, where we assume a "gluon mass" A/\/2, 
and we have written the form of the effective Lagrangian at q 2 ~ A 2 , integrating out 
the massive gluon, and truncating on dim = 6 operators. 

Upon Fierz-rearrangment of the interaction Lagrangian, keeping only leading 
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terms in 1/N C , eq.(77) takes the form: 

— 9 2 — — — — 

~ ^7^75^7^75^) (78) 

where tpL — (1 — r y 5 )ip/2, tp R = (1 + 75)^/2. Here r A are Pauli matrices acting upon 
the isospin indices. 

For the present analysis we will truncate eq.(78) on the pure Nambu-Jona-Lasinio 
terms, since the (vector) 2 and (axial- vector) 2 terms play no significant role in the chiral 
dynamics (they are associated with the formation of virtual p and A\ vector mesons 
in the model). Hence we take: 




We can solve the light-quark dynamics in large-iV in the usual way by writing an 
equivalent effective Lagrangian of the form: 

C = $W -m^-# L S^-^SV L -iA 2 Tr(StS) (80) 

where: 

S = i(7/ 2 +m a y (81) 

is an auxiliary field. We emphasize that at this stage E is a 2 x 2 complex field, so 
both a and 7r a are complex, (otherwise, with and ir real there would be unwanted 
contributions from (TE E) = E^ 7^ in integrating out E). 

Thus there is parity doubling at this stage, Im(a) is the fourth Goldstone boson 
associated with the U(l) problem, and Im(TT a ) is the + isotriplet. The restriction to 
real 7r a will emerge dynamically at very low energies, since the induced T^E^EE^E) 
term will lift the degeneracy of the Re(7r) and 7m(7r). We ultimately must add a 
det(S) + h.c. term to get rid of the Im(a) mode. 
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We now integrate out the fermion fields on scales A 2 > q 2 > fi 2 , keeping only the 
leading large- Nc fermion loop contributions. We use the massless fermion propagator, 
treating E as a classical background field. Thus we arrive at an effective field theory 
at the scale ji: 



C = 4>{i$ - m q )ip - g^L^tpR - gip^i/jL 
+Z 2 Tr(d l 0& t Yi) -Vp) 



(82) 



where: 



Z 2 = 

vp) = 



9 2 N 
16tt 2 



ln(AV^) 



I A 2 -^ 
2 8tt 2 



(A 2 - „ 2 ) 



Tr(£t£) - f^(A 2 - fi 2 ) Tr(m 9 E + h.c.) 



9 4 N 
16n 2 



ln(A 2 //i 2 ) 



(83) 



We see that Z 2 — > as /i — > A, reflecting the compositeness of the £ field. Let us 
now renormalize the £ field: 

£ -> (84) 

and we have the properly normalized effective Lagrangian at the scale \x (this is proper 
normalization for real a and 7r): 



- m q )4> - gip L T,ip R - gi) R Y)i\) L 
+ Tr(9 M S f ^S) -Vp) 



(85) 



where: 



Vp) 



m 2 Tr(S f S) - uo Tr(m g S + /i.e.) 
+ATr(S t SS t S) 



m„ = 



1 



1a 2 

2 



9 2 N f ,o 2 n 

8^ (A - " » 
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167T 2 = _ 2 

Aln(A 2 //z 2 9 

The effective Lagrangian is seen to be a linear a-model at scales /i < A. As the scale 
/t^Owe see that the theory is trivial, since g — > 0. However, these evolution results 
apply only to a scale /i corresponding to a mass scale for the fermion. Nonzero m q 
will block the evolution into the far infrared, but we will neglect this presently. The 
theory will develop a chiral instability (a constituent quark mass) provided that m 2 
becomes tachyonic (negative) at some scale /xq- By tuning the bare coupling constant 
g 2 we can put the model in a symmetric phase, m 2 > — > g 2 N /An 2 < 1, or in a chiral 
symmetry breaking phase: m 2 < — > g 2 N/4ir 2 > 1, where the critical bare coupling 
corresponds to m 2 = as /Jo — > 0. 

In the broken phase (ignoring m g ) the a field develops a vacuum expectation value 
(a) = f n = V2\m a \/ VX. We see that the renormalized a field develops a vacuum 
expectation value given by: 

«:-*(«^)(£-0-(f)(SH (87) 

In the broken phase we can then write a = f n + a, and the physical mass 2 of the a 
is readily seen to be m 2 = 2|m 2 |, while the fermion mass becomes m = \ j^g. Thus, 
using eqs.(86) to relate g 2 = A, we obtain the usual Nambu-Jona-Lasinio result: 
mj = 2m . 

The solution to the theory can thus be written as a chiral quark model in which 
we have both constituent quarks described by ip and the mesons decribed by S. In 
the broken phase it is useful to pass to a nonlinear cr-model and write: 

E^^exp^TrVVA) (88) 

and: 

C = -m q )ip-m ip L exp(in a T a /f w )^ R -m 'ip R exp(-iiT a T a /f 7T )tl) L 

+ Tr^E^E) + to Tr(m ? E + h.c) (89) 



-37- FERMILAB-Pub-93/059-T 



where m = \gfn is the constituent quark mass. Note, in our present normalization 
conventions that f w = 93 MeV. By a chiral redefinition of the fields, ip R — > £ip R 
and ip L — * ^tpL we arrive at the Georgi-Manohar Lagrangian (their eq.(2.9)) with 
qa = 1-0 (note that they fit Ga/Gv = (5/3)gA and obtain qa = 0.75, consistent with 
our large-N approximation). 

When the a and rr fields are slowly varying in space, the light quark propagator 
of the chiral quark model is given by (in terms of the unrenormalized fields): 

1 , 1 



^ - m q - g-a - ig-^n ■ r^j 

f j) + Tn q + gS 5 \ 



p 2 -Q 2 



(90) 



where we define: 



fr = (m q + #a/ 2 ) + #V/4; vr" = 7r a 7r a (91) 

s5= ^ + 4 757r ' r (92) 

In the broken phase we replace a = f T and E — > \f n ex. i p(m a T a / f n ). For future ease 
of writing we can often replace ga/2 = ga/2 + m q ^fZ^ since it easy to restore the 
explicit chiral symmetry breaking quark mass terms. 

(ii) Schematic Discussion of a Parity Doubled Nucleon 

Consider a "nucleon" doublet iV with the SU(2) L x SU(2) R assignments N L ~ 
(|,0), N R ~ (0, |). Also, we introduce a partner, K, of opposite parity with assign- 
ments K L ~ (0, |), ~ A typical renormalizeable linear a-model effective 
matter Lagrangian (not including the E kinetic and potential terms) is then: 

C = ~Ni(J) N + ~Ki(/) K 

-MJTlENr - M 2 K L tfK R - M N L K R - M' N R K L + h.c. (93) 

Parity symmmetry requires M = Mq. We consider the special case Mi = M 2 = M, 
which is the analogue of our model, but this is not generally required by symmetries. 
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Now perform the redefinitions, N L -> £N L , K L -> ft# L> iV R -> ^Nr, K r -> £,Kr. 
Thus, the Lagrangian becomes: 

-MNN - MKK - .U„A7\ - M ~KN + /i.e. (94) 

Upon diagonalizing, the mass eigenfields are just (N±K)/y/2, with mass eigenvalues 
M ± M . We can decouple the heavier state by setting (N + K)/ \/2 = 0, whence the 
light effective Lagrangian for Q = (N — K)j \/2 is: 

£ = Q{i$ + V )Q - (M - M )QQ (95) 

We see that qa = 0. Hence, g A is not generally of order unity as is the case of a 
non-parity doubled nucleon. (this is also a consequence of the special case Mi = M 2 ; 
more generally qa = sin(26 ) ) where 9 is the mass mixing angle). With g — the 
only nontrivial Goldberger-Treiman relation refers to the pionic transition amplitude 
between the ground state, Q, and the parity partner. 
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